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Abstract. We calculate the optically thin synchrotron emission of fast electrons and positrons in a spiral stationary 
magnetic field and a radial electric field of a rotating relativistic strongly magnetized force-free jet consisting of 
electron-positron pair plasma. The magnetic field has a helical structure with a uniform axial component and a 
toroidal component that is maximal inside the jet and decreasing to zero towards the boundary of the jet. Doppler 
boosting and swing of the polarization angle of synchrotron emission due to the relativistic motion of the emitting 
volume are calculated. The distribution of the plasma velocity in the jet is consistent with the electromagnetic field 
structure. Two spatial distributions of fast particles are considered: uniform, and concentrated in the vicinity of 
the Alfven resonance surface. The latter distribution corresponds to the regular acceleration by an electromagnetic 
wave in the vicinity of its Alfven resonance surface inside the jet. The polarization properties of the radiation have 
been obtained and compared with the existing VLBI polarization measurements of parsec-scale jets in BL Lac 
sources and quasars. Our results give a natural explanation of the observed bimodality in the alignment between 
the electric field vector of the polarized radiation and the projection of the jet axis on the plane of the sky. 
We interpret the motion of bright knots as a phase velocity of standing spiral eigenmodes of electromagnetic 
perturbations in a cylindrical jet. The degree of polarization and the velocity of the observed proper motion of 
bright knots depend upon the angular rotational velocity of the jet. The observed polarizations and velocities of 
knots indicate that the magnetic field lines are bent in the direction opposite to the direction of the jet rotation. 

Key words, radiation mechanisms: non-thermal - magnetic fields - galaxies:jets - BL Lacertae objects: general - 
quasars: general 



^ 1. Introduction 

5-j \ 

The jets from Active Galactic Nuclei are not uniform. They consist of a number of bright knots with fainter emission 
in between. It is well established that the nature of the radio emission from jets on kiloparsec and parsec scales is 
synchrotron radiation of relativistic nonthermal particles (Begelman et al. 1984). There is mounting observational 
evidence that jets in many AGNs are electron-positron rather than electron-proton. The spectral model of the jet 
radiation in a broad frequency range (from X-rays to radio), is described by Morrison & Sadun (1992), who pointed out 
that the 3C273 jet must be dominated by electron-positron pairs. Reynolds et al. (1996) used early VLBI observations 
of the M87 jet and self-absorbed core. They used the theory of synchrotron self-absorption to put a constraint on the 
electron density in the core and concluded that the core of M87 is likely to be e + e~ pair plasma. A similar method was 
used by Hirotani et al. (1999) and Hirotani et al. (2000) to reveal that the components in the jets of 3C279 and 3C345 
are also dominated by pair plasma. Observations of circular polarization of radio emission from 3C279, 3C273, 3C84 
and PKS0528+134 (Wardle et al., 1998) support the same conclusion. However, there are also indications that jets in 
Optically Violently Variable quasars cannot consist solely of e + e~ pairs because they would produce much larger soft 
X-ray luminosities than are observed. On the other hand, models with jets consisting solely of proton-electron plasma 
are excluded, since they predict much weaker nonthermal X-radiation than observed in Optically Violently Variable 
quasars (Sikora & Madejski, 2000). The in situ acceleration of relativistic particles seems unavoidable because high 
synchrotron losses prevent relativistic particles from streaming out of the central engine to the distances of bright knots 
(Begelman et al. 1984; Roser & Meisenheimer 1991). There are several possible mechanisms of particle acceleration 
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inside knots. The acceleration can be by shocks (Blandford & Eichler 1987), by turbulent electromagnetic fields 
(Manolakou, Anastasiadis & Vlahos 1999) and by magnetic reconnection (Romanova & Lovelace 1992; Blackman 
1996). Another possible mechanism of particle acceleration is the acceleration by organized electromagnetic fields 
which are the eigenmodes of the jet as a whole. These modes can be excited as a result of non-stationary processes 
in the central engine. The energy can be transmitted by such global waves along the jet to large distances without 
noticeable attenuation. 

The acceleration is most effective if resonances of the waves with the mean stream take place (Beresnyak et al. 
2002). The resonance with Alfven waves occurs when u>' = k'^c, where lo' and k'^ are the frequency and the projection 
of the wave vector on the magnetic field in the plasma rest frame. When the specific energy density and pressure of the 
plasma are much less than the energy density of the magnetic field, The Alfven velocity is equal to the speed of light 
in vacuum c. The resonance condition is fulfilled on a specific magnetic surface and, in the case of a cylindrical jet, 
at a definite distance from the axis. In the vicinity of that surface the magnetic and electric fields of the wave reach 
large magnitudes. The resonances occur for nonaxisymmetric modes of global waves and the position and width of 
the resonances depend on the structure of the ordered helical magnetic field. Both electrons and positrons are subject 
to drifting out of the region of the strong electric field of the wave near the resonance due to the non- uniformity of 
the electric field. When moving to a region of weak electric field, particles gain energy. This energy is taken from the 
energy of the wave, which is an eigenmode in the jet, and the eigenmode decays. Beresnyak at al. (2002) considered 
this mechanism of particle acceleration. Drift approximation was used to describe the motion of a particle in the 
electric field of the wave. Synchrotron losses and isotropization of the particle distribution by plasma instabilities were 
taken into account. Plasma instabilities are excited by accelerated particles if the distribution function is anisotropic in 
momentum space, and act to conceal any anisotropy in the particle distribution. Acceleration process and synchrotron 
losses taken together yield a power law energy spectrum of ultra-relativistic electrons and positrons with an index 
varying between 2 and 3 depending upon the initial energy of the injected particles. This is consistent with the typical 
spectral indices of radio emission observed in parsec-scale jets. 

The polarization of synchrotron emission in knots is very sensitive to the geometry of the helical large scale magnetic 
field. Indeed, VLBI observations provide evidence for the polarized emission in knots and in the space between knots 
(Gabuzda 1999a; Hutchison et al. 2001). This is indicative of the existence of a large scale magnetic field all over the 
jet rather than concentrated in separate knots. In the general case a helical magnetic field has poloidal and toroidal 
components. The toroidal component arises due to the axial electric current in the jet. If the current inside the jet 
is closed (no net jet electric current), that is the electric current flows in different directions at different radii in the 
jet, then the toroidal magnetic field must be small at the boundary of the jet and only the axial component of the 
magnetic field will remain. The change of the direction of the magnetic field from toroidal close to the jet axis to axial 
on the periphery was inferred from the VLBI polarization measurements in blazar 1055+018 (Attridge et al. 1999). 
The concept of an intrinsic large-scale helical magnetic field is the simplest model of the jet and does not require 
additional processes of amplification or generation of azimuthal and axial components of the magnetic field in shocks 
or shearing layers at the jet boundary. 

The polarization of synchrotron emission of parsec scale jets is being studied by few groups (Gabuzda 1999b, 
2000; Leppanen at al. 1995; Lister 2001). We use observational data to compare with our calculations of polarization 
in a model of a force-free strongly magnetized relativistic rotating jet. We calculate the polarization of synchrotron 
radiation in the helical field geometry (Sect. [21 • In Sect. we expand our calculations to include the Doppler boosting 
effect and swing of the polarization angle due to relativistic motion and rotation of the plasma in the jet. We also 
compare the observed apparent velocities of knots with the velocities given by our model of propagation of spiral 
perturbations in a cylindrical jet. As it was shown by Istomin & Pariev (1994, 1996) there exist so-called standing 
modes among eigenmodes in jets (w gr0 up = 0) which do not propagate with finite velocity along the jet but arc only 
subject to diffuse spreading and, therefore, their amplitude is maximal. The phase velocity of these modes is greater 
than the speed of light. The crests of the wave move along the jet with super luminal velocity causing acceleration of 
the particles on the Alfven resonance surface. These regions can be the objects which are seen as bright knots with 
typical sizes of the order of the wavelength of the standing wave, which is about the radius of the jet. 



2. Synchrotron emission. 

Let us now consider synchrotron emission of accelerated particles in a cylindrical force-free jet with a spiral magnetic 
field. For cylindrical jet the stationary configuration of fields is (Istomin & Pariev 1994): 

B = B z e z + B e ; = ± B z , (1) 
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where r, 4>, z are cylindrical coordinates. The profile of the poloidal magnetic field B z is taken to be uniform across 
the jet. The velocity of plasma can be written as 

u = KB + ft F re , (2) 

where fi F is the angular rotational velocity of the magnetic field lines and KB is the component of the velocity parallel 
to the magnetic field. The electric field is induced due to the motion of highly conducting plasma across the magnetic 
field and is equal to 

E = r [e^xB]. (3) 

In components 

B = (Q,±n F r/c,l)B z , (4) 
E= {-n F r/c 1 7 0)B Z7 (5) 
u= (0,n F r(l±KB z /c),KB z ). (6) 

Here £l F and K are functions of r, and B z does not depend on r. The stationary electric field, which has an absolute 
value of fl F rB z /c, is not small compared to the magnetic field when Q, F r ~ c. The rotation velocity Vi F r can exceed 
the speed of light; nevertheless u remains less than c due to the existence of the predominantly toroidal magnetic 
field. This can be achieved by appropriate choice of the parameter K . The natural requirement that the total poloidal 
current is closed inside the jet and that the total electric charge of the jet is zero implies that fl F must vanish at the 
boundary of the jet r — R (Istomin & Pariev 1994). Besides this requirement, the dependence fl F — il F (r) should be 
determined by the conditions at the jet origin. An arbitrary parameter K(r) in expression for the plasma velocity 
determines the radial profile of the longitudinal velocity. K(r)B z is the velocity along the magnetic field which is not 
related to the rotation. We choose a velocity u which minimizes the kinetic energy of the plasma in the stationary 
reference frame 

(0,^r, T (^V) 2 / C ) 

l + (fl F r/c) 2 ' [) 

It is probable that the plasma moves with that velocity in real jets originating near black holes (Frolov & Novikov 
1998). The geometry of the magnetic and electric fields and the flow in the jet are illustrated in Fig. ^ 

In accordance with our theory of particle acceleration inside a force-free relativistic jet (Beresnyak et al. 2002) we 
consider the distribution function of emitting particles to be isotropic in momentum and power law in energy 

dn = K e £- c d£dVdn p (8) 

Here dn is the number of particles in the energy interval £,£ + d£, dV is the elementary volume, gK7 p is the elementary 
solid angle in the direction of the particle momentum p, K e — K e (r), £ = constant. We will assume that ultrarelativistic 
particles are distributed according to Eq. (|SJ) in the frame moving with the plasma with the velocity u, Eq. JJJ. One 
needs to know the distribution of the density of relativistic particles over the jet radius K e (r). Several mechanisms for 
particle acceleration inside jets have been considered previously. The possible processes are acceleration by the shocks 
in the supersonic flow of mass-loaded jets (Blandford & Eichler 1987) and acceleration due to magnetic reconnection 
near the neutral layers of the poloidal magnetic field (Romanova & Lovelace 1992). Usually, the distribution of fast 
particles inside the jet was assumed to be uniform (Laing 1981) or smoothly varying with the radius in the jet (Konigl & 
Choudhuri 1985). According to our mechanism of particle acceleration (Beresnyak et al. 2002), relativistic particles 
are accelerated close to the Alfven resonance surface r — r a and are concentrated near the same surface. The position 
of the Alfven resonance in a force-free jet is determined by the condition 

kc~LU + 2mQ F (r A )- ^^a)^ ( w + fcc ) = 0; (g) 

where uj = uj m (k) is the eigenfrequency of the perturbation of the magnetic field cx exp(— iujt + ikz + imcj)) with the 
wavenumber k along the axis of the jet and azimuthal number m. We will consider two cases for the spatial distribution 
of particles, namely, uniform in the jet volume and concentrated close to the Alfven resonant surface r — Ta- 

The influence of the plasma on the emission of relativistic e~ and e + is neglected in the present work. Absorption 
and reabsorption are also not taken into account. Electrons and positrons are accelerated in the same way near 
Alfven resonance and have the same distribution function. Therefore, the Faraday rotation effect is absent as well as 
intrinsic depolarization inside the jet. VLBI polarization measurements are done at frequencies in the range 1 GHz 
to 100 GHz. Synchrotron self-absorption is very small at GHz frequencies and higher frequencies. For B ~ 10~ 2 G, 
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n ~ 0.1cm 3 , and C, ~ 2.5 an order of magnitude estimate for the synchrotron self-absorption coefficient is k s — 
(5 • 10 3 pc) _1 [y/l GHz)~ 2 ~''/ 2 (e.g., Zhelezniakov 1996). The corresponding optical depth through the 1 pc jet is only 
2 ■ 1CP 4 . The typical energy of particles emitting in this frequency range in 1CP 2 G magnetic field is 100 MeV to 
1 GeV. Conversion of linear polarized radio waves into circular polarized (Cotton-Mouton effect) is also small. For 
both thermal and relativistic power law particle distributions with £ w 2.5 and a minimum Lorentz factor ~ 10 the 
estimate for the conversion coefficient is k c — (500 pc) -1 {v/1 GHz) -3 (Sazonov 1969). This is too small to influence 
the transfer of linear polarization inside the jet, but can be the source for a small, at the level of a fraction of per cent, 
circular polarization. We focus on linear polarization in this work. 

At first, we neglect Doppler boosting, i.e. assume that u = for the purpose of the calculation of the synchrotron 
radiation. However, the magnetic field will still be given by formulae with nonvanishing f2 F . In Sect.OHwe take into 
account Doppler boosting with velocity u consistent with the field configuration according to expressions and 
see that the proper consideration of the Doppler effect leads to a marked change in polarization. When no Doppler 
effect is considered, the Stokes parameters are obtained by integration along the line of sight passing through the 
region filled with relativistic particles inside the jet (e.g., Ginzburg 1989). Let us denote the angle between the axis of 
the jet (the direction of the axial magnetic field B z ) and the observer by 9, such that 9 = when B z points exactly at 
the observer, 9 = tt/2 when the line of sight of the observer is perpendicular to the direction of the jet axis, and 9 = tt 
when the B z points exactly away from the observer. The azimuthal angle <j> is measured from the plane containing the 
observer and the jet axis; <f> increases from to 2tt in anticlockwise direction if viewed from the positive direction of 
axial magnetic field B z . Let us denote the distance between the line of sight and the jet axis by h and the radius of 
the jet by R. We will take h as a signed quantity: h > for lines of sight passing in the sector < 4> < tt, h < for the 
lines of sight passing in the sector tt < <f> < 2tt. We will parametrize integration along the line of sight by integration 
in </> keeping the value of h constant. We define Stokes parameters in the usual way with reference to the rectangular 
coordinates in the plane of the sky with the x-axis oriented along the projection of the jet axis on the plane of the sky 
and the y-axis orthogonal to the x-axis in the anticlockwise direction. This geometry and the angles are illustrated 
in Fig. ^ Considering the geometry described above we obtain the following expressions for the Stokes parameters 

02 

1 = C ~rrr k ^ I l^l (c+1)/2 * 2 . ^ (w) 

C + 1 J sinfl sin <p 

01 

r, \ /"V±| (C+1)/2 [B 2 COS 2 <j>~(B z Sm9 + B^Sm ( f ) CO S 9) 2 }h 

Q = k(v) r^-ja . - . 2 7 #> ( n ) 

J \B±\ sin # sin <p 



01 



U = -2k{v) / 1 j ' — B^ cos (f)(B z sin 9 + B sin </> cos 9) ^— d<f>, (12) 

' l-DJ-l sin 6> sin A 



V = 0, (13) 

where Bj_ is the component of the magnetic field perpendicular to the line of sight. The integration limits are <j>\ — 

7rsgn/i — arcsin — , 2 = arcsin — . The value of k{v) is a function of the observed frequency v of the radiation 
R R 

(K- 1\ r, /3C + 7\ e 3 r 3e 1 (C-1)/2 _ lt ._^ ln K, 



4 V 12 / V 12 



27rm 3 c 5 



-(C-l)/2fl£ 



where D is the distance between the source and the observer, e and m e are the charge and mass of an electron, and 
r is the Euler gamma-function. 

We denote the position angle of the electric field vector in the plane of the sky by \- The angle x is measured 
clockwise from the direction parallel to the projection of the jet axis on the plane of the sky and < x < t ( see Fig.^). 
For ultra-relativistic particles 1^ = 0, and the emission has only linear polarization. The degree of polarization of the 
observed radiation is expressed as n = \JQ 2 + U 2 /I, the resultant position angle of the electric field Xrcs measured 
by the observer is found from 

Q U 

COS2Xrcs= 7= tt9 , Sin2Xrcs = —====, 0<Xrcs<7T. 
y/Q 2 + U 2 y/Q 2 + U 2 

It is seen that when one replaces the variable 4> by tt — tj> in expression (|12fl . U changes its sign and, therefore, U = 0. 
Consequently, if Q > then Xxes — 0, if Q < then Xrcs = tt/2. Thus, the observed electric field vector is either 
parallel or perpendicular to the jet axis on the plane of sky. This may account for the bimodal distribution of the 
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observed x rcs (Gabuzda et al. 1994). Now let us integrate expressions (|10fl and for I and Q over parameter h in 
order to obtain the total radiation flux and polarization per unit length of the jet image on the sky. Assuming uniform 
emissivity over the jet volume, B z = constant and — B z il F r/c expressions for / and Q from the unit jet length 



become 

<t>2 



/=^f%) fdhf \B ± \W* , k d0, (14) 
C + 1 J J sm sin <p 

—R <Pi 

Q = k{ V ) jdh J\Bj<-M ^ C ° S2 = iC f n9 + " f r BMcoBfl) W # , (15) 
J J c 2 sin 8 sin </> 

-R <t>i 

For C = 3 expressions (11411 and ((15(1 are reduced to the form 



R 





R 



n p2 r 2 ( 1 

sin 2 6 H = — 1 - - sin 2 i 



// O-F^ 2 \ 
— ^-2 rdr. (17) 

o ^ ' 

From the above we see that in the case £ = 3 the sign of Q does not depend on 8 since sin# > for any < 8 < it. 
The sign of Q is determined by the function fl F (r) only. 

Let us choose a particular dependence of Q (r) satisfying the requirements f2' F (i?) = and dft F jdr — at r = 0. 
At present, there is no well established theory which would predict £l F (r). Therefore, we choose the simplest monotonic 
profile satisfying the physical requirements above: 

"'-"sHs)")- (18 > 

where f2 is a dimensionless parameter. Although the exact position of ta, the ratio of B^jB z at ta and the velocity of 
the flow in the vicinity of ta (properties which determine the observed polarization) are all dependent on the particular 
choice of ft F (r) profile, it is the magnitude of the angular rotational velocity which influences all these parameters 
most significantly. A qualitative difference may appear if the profile of £l F is non-monotonic, such that many Alfven 
resonant surfaces will exist inside the jet. In this work we limit our consideration to variation of parameter f2 only, 
but keeping the radial profile of tt F monotonic and given by expression (|18|l . Then, for £ = 3 using expressions l|16[) 
and HI 7(1 one obtains the degree of polarization II 

(— - 1^ sin 2 ' 



dr, (16) 



4 • 2n " 2 (■> 1 • 2 

1~2 ~ ' 



If |fi| > 2-^/6, then Q > and Xrcs = 0, i.e. the electric vector of the polarized emission is parallel to the projection 
of the jet axis on the plane of the sky. This does not mean, however, that the orientation of the magnetic field in 
the jet is purely perpendicular to the jet axis. Clearly, one can see from Eq. Q that the poloidal magnetic field is 
comparable to the toroidal field. If < < 2^/6 then Q < and Xres = tt/2, i.e. the electric vector of the polarized 
emission is perpendicular to the projection of the jet axis on the plane of sky. Again, this does not mean that the 
orientation of the magnetic field in the jet is purely parallel to the jet axis. The toroidal magnetic field vanishes only 
if Q = 0, otherwise the toroidal magnetic field is comparable to the poloidal magnetic field. This simplified analytical 
calculation shows that one needs to be very careful when interpreting the polarization measurements in terms of the 
orientation of the magnetic field. 

Eqs. ((14(1 and l(15|) were integrated numerically for the dependence tt F (r) given by expression l(18|) and £ = 2.1, 
which corresponds to the power law index in the spectrum of synchrotron emission of (£ — l)/2 = 0.55, the mean 
observed value for extragalactic jets (Scheuer 1984). The results are presented in Fig. El Positive values of the degree 
of polarization correspond to Xrcs = 0, negative to % rcs = 7t/2. In the case II = the radiation is entirely unpolarized. 
Note also that in the nonrelativistic consideration 11(180° — ff) = 11(6') for any value of 8. 

Istomin & Pariev (1996) presented a hypothesis that knots observed in jets can be the manifestation of the "standing 
wave" phenomenon. If one knows the dispersion curves uj = uj{k) and wavenumbers fc m ; n at which duj/dk = for 
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Table 1. Polarization degrees and velocities of knots in BL Lac objects and quasars 



Class 


Source 


n,% 


v v i s h/c 


References 


BL Lac 


0138-097 


8.8 


8.5 


1 


BL Lac 


0454+844 


8.5 


0.92 


1 


BL Lac 


0735+178 


-14. 


8.5 


1 


BL Lac 


0828+493 


-12.5 


9.7 


1 


BL Lac 


0851+202 


-23.3 


4.3 


1 


BL Lac 


0954+658 


17.2 


8.8 


1 


BL Lac 


1418+546 


5.6 


3.2 


1 


BL Lac 


1732+389 


-25.9 


17.8 


1 


BL Lac 


1803+784 


27.5 


3.2 


2,3 


BL Lac 


1823+568 


14.9 


6.2 


1 


BL Lac 


2007+777 


5.6 


4.5 


1 


BL Lac 


2200+420 


11. 


5.5 


1 


quasar 


0212+735 


1.4 


8.6 


2,3 


quasar 


0836+710 


2.2 


21.7 


2,3 


quasar 


0906+430 


7.5 


5.3 


2,3 


quasar 


1055+018 


9.5 


9.23 


4 


quasar 


1308+326 


3.8 


16.8 


5 


quasar 


3C345 


11. 


29.3 


2,3 


quasar 


1642+690 


13. 


14.0 


2,3 


quasar 


1928+738 


9.3 


11.5 


2,3 


quasar 


3C380 


24.4 


10.6 


2,3 



RefeRences.— (1) Gabuzda et al. (2000), (2) Lister et al. (2001), (3) Lister (2001), (4) Lister et al. (1998), (5) Gabuzda & 
Cawthorne (1993) 



perturbations, one can calculate the velocity of moving crests of the "standing wave" as u p has = ^min(^min)/&min 
and corresponding observed velocity v v i s — phas — . Note also that changing the sign of Bj, in Eq. ffl is 

1 - fphasCOSfcyc 

equivalent to changing the sign of f2 F , since / and Q depend only on B z and B^ and not on u or E r . In Fig. [3 we 
present the relations of the polarization Id of the jet emission and v V is obtained in the case of O changing from to 
100 for some angles 9. We also plotted observational points for jets in BL Lac objects and quasars for which we were 
able to find both polarization and proper motion measurements in the literature and which are listed in Table 1. For 
each object we averaged the degree of polarization and observed velocity over all bright knots observed in a particular 
object excluding a few knots, for which peculiarity was indicated in the original references. This procedure might be 
very questionable although because of the lack of observations of the emission from the space between knots on parsec 
scale we had nothing better. The polarization measurements also depend on the frequency of observations. Specifically, 
different observations of quasars give very different values of the polarization of the knots at centimetre and millimetre 
wavelengths. Numbers for v v - ls were calculated assuming a Hubble constant H = 65/ikm/s • Mpc 

In the case of accelerated particles concentrated very close to the Alfven surface r = ta expressions (|14|) and 
1)15(1 can be reduced to integrals in (p only, while r = ta- Since for axisymmetric modes (m — 0) always |w m in| > 
|&min|c (Istomin &; Pariev 1994) one can see from the expression © that the Alfven resonance does not exist for 
axisymmetric modes. Therefore, there should be no particle acceleration for this mode. The next most important 
modes are m = 1 and m = — 1. Since none of the quantities fc m i n , Wmm or II changes when m changes to —m and 
Q F to — f2 F , it is enough to consider only positive m = 1 and all values of £l F (positive and negative). Polarization 
Id versus f p has for this case is plotted in Fig. as well as observational points. We plotted i> p has only for the first 
non-axisymmetric fundamental mode m = 1 which is expected to be excited more easily than higher modes with 
m > 1. Observations of the parsec scale jet in BL Lac object 0820+225 (Gabuzda et al. 2001) show the S-like shape 
which can be attributed to the mode m = 1. The polarization curves are plotted only for values of Q at which the 
Alfven resonance r = ta is present inside the jet. As previously, Q F (r) is given by expression (|18fl and C — 2.1. It is 
seen that in the case of uniform particle distribution across the jet (Fig.EJ), curves for changing from 120° to 70° are 
located in the region where points from observations of BL Lac object (crosses) and quasars (squares) are. However, 
the polarization curves for emitting particles concentrated on the Alfven resonance surface (Fig.^J, as is the case for 
our acceleration mechanism, do not match the observations of quasars. It is seen from Fig. [3] that the electric field 
vector is always oriented perpendicular to the jet axis unlike the case with uniform distribution of relativistic particles 
across the jet. The reason for this difference is that the Alfven resonance surface exists only for moderate values of 17, 
when the toroidal magnetic field is not strong enough to cause the x res to become oriented parallel to the projection 
of the jet axis on the plane of the sky. 
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3. Relativistic Doppler boosting effect on polarization 

Now let us make more realistic calculations taking into account Doppler boosting and the presence of a strong electric 
field. As was first pointed out by Blandford & Konigl (1979) the effect of relativistic aberration of light causes a swing 
of the direction of polarization of synchrotron emission when the source is Doppler boosted. This is due to the fact 
that the electric vector of the polarized component of the wave must remain orthogonal to the wave vector. We will 
denote values measured in the plasma rest frame by primes. Let us find the transformation rules for the polarization 
tensor J a p = < E a E*j > from the local reference frame of the small emitting volume moving with velocity u to the 
observer reference frame. Here < . . . > is the averaging over realizations of electric wave fields. The Stokes parameters 

are components of the polarization tensor J a p = 1^ Q^) ^ G denote the local Cartesian coordinates in 

the source frame by x'y'z' and the Cartesian coordinates in the observer frame by xyz. The transformation can be 
done by directly using the transformation formulae for the electric and magnetic fields of the wave (Cocke & Holm 
1972). The result is 



l-u 2 /c 2 

Ja(3 = J a '8'Ti / rv>- (^) 

P (1 - u/ccosS) 2 



Here we denote by 5 the angle between velocity u of the source and the direction toward observer k measured in the 
observer frame; correspondingly, 5' is the same angle but measured in the source frame. Thus, the polarization ellipse 
in the observer frame is obtained from the polarization ellipse in the source frame by rotation by the angle 5 — 5' 
in the plane containing the direction to the observer k and u. Therefore, the degree of polarization is not changed, 
which is obvious. The Stokes parameters for the radiation of one particle are expressed through the radiation flux 
with two main directions of polarization p^) and p^) by the formulae I' = p^} + p$ , Q' = {p^) —p^))cos2yf , 
U' = {p^) — p^)sm2%', where x' is the position angle of the electric vector of polarized radiation in the reference 

frame comoving with the emitting volume (Ginzburg 1989). Then, p^) and p^} are transformed to the observer frame 
as 

„(1) = J^— 1 - U2 /° 2 D (2) = J 2 ) 1 - u2 /° 2 

Pv Pv ' dv {1-u/ccosS) 2 ' Pv Pv ' dv (l-w/ccos<5) 2 ' 
Therefore the intensity in the observer frame integrated along the line of sight is 

/(i/, k) = / (pW ry, e', l, x ', vo + p {2) K e', l, x', /)) x 



dv' 1-uVc 2 ,- C dV , 

~TT\ / Tvi Ke£ -fnrdS dflp.dL, 

dv (1 — u/ccosdy dv 

where dL is the element of length along the line of sight, \' is the angle between the momentum of the particle and the 
magnetic field vector, ip' is the angle between the direction of radiation and the cone formed by the velocity vector of 
the particle rotating around the magnetic field line, both measured in the frame comoving with the fluid element. We 
will integrate over the momenta in the plasma rest frame while integrating over the space coordinates in the observer 
frame. Integration over the solid angle f2 p / in the momenta space in a comoving frame may be done analogous to the 
usual case of synchrotron emission of ultra-relativistic particles, when there is no electric field E = u x B. Integration 
over f2 p / and £' leads to (Ginzburg 1989) 

/ \ 

t, , n C + 7/3 , , . f I 1 — u/ccos5 \ . c+i dL 

I(u,k) = ^—^-k(v) [^=J==) (B'sinx')^- j— -=. (20) 

C+l J y^/l-u 2 /c 2 J l -u/ccosd 

1 — u / c COS S 

Here v' = v — =■, B' is the magnetic field in the emitter frame. The expressions for Q and U can be obtained 

yj\ — u 2 /c 2 

in a similar way. The result is 

\ 

, , 1 — u/ccos5\ . ,. c+i cos2ydL 

= W] \ -JY=uW ) { ^ 2 1-ti/ccosJ ' (21) 

, ii l-u/ccos5\ . .. <+i sin2vdL 

U(v,k) = kM / \ ' (B sin y )~ r (22) 

1 ' V 'J \ ^l-u 2 /c 2 J K A7 I-m/ccos.5 v ' 
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The component of the magnetic field perpendicular to the momentum of the particle in the comoving frame, B' sinx', 
is obtained by applying the transformation rules 



B' 2 sin 2 X ' = B 2 +B 2 -E 2 - 

[— E r (c 2 — u 2 ){u z sin <p sin 9 + u,p cos 8) + u(c cos 8 — u)(it^£f^ + u z B z )} 2 

u 4 (c — ucosS) 2 



(23) 



The expression for the position angle of the electric vector of the polarized radiation from a jet element in the observer 
frame \ is obtained from yf by projecting the polarization ellipse in the plane perpendicular to k' on the plane 
perpendicular to k. After some geometrical considerations the result for x can be expressed as follows 

X = Xu + °, (24) 

where angle x.u can be found from 

E r (u — ccos S)(u z sin <j)sm8 + us cos 9) — u sin 2 8{usBs + u z B z ) 

tanx u = — : — — -r , (25) 

k, r u sm cos 0(c — ucosd) 

< Xu < 

and angle a, < a < 2it, can be found from 

us sin <p cos 6 + u z sin 9 Us cos tfi 

cos a = : — , sin er = —. — — , (26) 

u sm d u sm d 

where the expression for 6 is 

u cos (5 = — us sin 9 sin</> + u z cos 9, < 8 < it. (27) 

Note, that by simply putting u = in Eqs. I|20II27|I one recovers the formulae for the emission from a jet with B& — 
and E r = . In order to recover Eqs. Q14H15JI one needs to keep Bs and E r unchanged but make no distinction between 
the reference frame of the emitting element in the jet and the observer frame. For uniform particle distribution the 
intensities from a unit length of the jet are 

( + 7/3 k(v) [ R , f 27r , , /l-u/ccosA - * 3 " |S'sinx'| (c+1)/2 

1 = — ■ — - I r dr 



C + l sin^7 




Q = —» f Tdr ** V H'"7il ) |£'sin X 'l^/* — «2S, 
sm 9 J Q J \ V 1 _ u l c J I -u/ ccos 5 

C-l 

u= m f R rdr /%( w£^£) 2 |^ sinx '|(c + D/ 2 . sin2 x 



o Jo \ x/ 1 - u2 / c2 / l-u/ccos<5 

When the particles are located close to the Alfven resonance surface in the region \r — Ta\ < A and A <C R, the 
integration in formulae (|28|l reduces to integration in <f> only while r = Ta- It can be checked that under the change 
of (j> to —cj) in Eqs. I|28|l the value of Q is not changed, and the sign of U is reversed. Thus, as a result of integration, 
U = and either Xrcs = or % rcs = tt/2. This result of the bimodality of the values of Xres is the same as in Sect. |2J 
where we do not take into account Doppler boosting and electric field. 

We perform numerical integration of Eqs. I|28l) using formulae l|23|l - H27|) . As in Sect. we take B z — constant, 
Bs = ±£l F rB z /c, and take the dependence fl F (r) given by expression (|18(l . 

The results for uniform particle distribution do not match the observations. Taking into account the presence of 
the electric field and Doppler boosting changes the polarization of the synchrotron emission essentially. The apparent 
direction of the magnetic field, as it is usually derived from polarization measurements (just perpendicular to the 
electric vector in the polarized component of the radiation), is now transversal to the jet axis only for small angle 9 
and large fl. For most inclination angles 8 and angular velocities Q, the apparent magnetic field is longitudinal (Fig.QJ. 
One can also see from Fig. where we plot IT versus v v - ls for Bs = +£l F rB z /c and for the particles concentrated 
near Alfven resonance surface, that our simulation does not match the observational data from Table 1. However, for 
Bs = —Vl F rB z /c we are able to reproduce the observations well enough, if the particles are distributed close to the 
Alfven surface and m = 1 orm = — 1. In both cases, with and without considering the Doppler boosting effect and 
the electric field, II is not changed when one changes the sign of fi, but II is changed when one puts Bs = —U F r/cB z 
instead of Bs = Q F r/cB z . Let us list all transformations which do not change the equations for the radial displacement 
£ r in the perturbations (Istomin & Pariev 1996) or lead to complex conjugated equations for £*: 
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1. 


LO — i 


• -w*, fc -> 


-jfc, n -> -0; 


2. 


LO — i 


■ -w*, fc -> 


— k, ni — > —rn; 


3. 


n - 


> —51, to — > 


— to; 


4. 


fc -> 


-fc, 




5. 


s,- 


' — -B z . 





-n F rB z /c; 

So, if the eigenvalue stability problem for the radial displacement has one set of parameters listed above as a solution, 
then the parameters transformed according to the rules 1-5 above will also be a solution of the eigenvalue problem for 
the radial displacement for the same £ r or for £*. Therefore, fc m i n and w m ; n can be transformed according to the rules 
1-5 above and the transformed fc m ; n and w m ; n will also correspond to the "standing wave" wavenumber and frequency. 
The value of r& is also not changed under transformation of the parameters listed above. The degree of polarization II 
is neither influenced by changing the sign of B z nor by changing sign of tt. The crucial fact is that a sign change in 
the relation — ±fl F rB z /c leads to a sign change of i> p has = <^imn/k m m, whilst all other four transformations listed 
above change neither II nor i> p has- Therefore, for the purpose of finding the correlation between II and v V is, which can 
be verified observationally, all different cases are reduced to considering two signs in the relation B^ = ±f!, F rB z /c, 
but we can limit ourself only to the consideration of B z > 0, m > and all O for which r& is located inside the 
jet. Using formulae H23 |l -I|27 |l and expressions JSJ for u z and u^, one can check that the Stokes parameters remain 
unchanged when one changes the sign of B^ and replaces angle 9 by 7r — 9. Also, a simultaneous sign change of w p has 
and change of 9 to 7r — 9 leads to v v is being transformed into — v v - ls . Therefore, the dependence of II on v V i S in the case 
£?0 = —Q F rB z /c can be obtained from the same dependence in the case S = +n F rB z /c by reversing the sign of 
^vis - * — Vvis and relabelling each curve from 9 to n — 9. 

We plotted the degree of linear polarization II versus u V is for to = 1 in Figs. 0] and for B^ — +fl F rB z /c and 
Bff, = —VL F rB z /c correspondingly. The Alfven resonance exists for 0.7 < Q < 4.9. Therefore, the value of belongs 
to this interval along the curves in Figs. 01 and [SJ 

We see that the curves in Fig. 01 do not match the observational data for knots in jets, while those in Fig. do. 
The dependence on the power-law index £ of emitting particles is weak. Our choice of ( = 2.1 and £ = 3 is motivated 
by the maximum and minimum values of £ obtained in the process of the formation of the spectrum of the particles 
accelerated near the Alfven resonance surface (Beresnyak et al. 2002). 

We can suggest the following explanation of the fact that knots are observed more often in jets with B^ = —Q F rB z /c 
than in jets having the opposite winding of the magnetic field. If we look at the models of the formation of magnetically 
driven jets from accretion discs (see, e.g. Pelletier & Pudritz 1992) we find that the magnetic field lines are always bent 
in the direction opposite to the direction of disc rotation as soon as the velocity of the outflow is directed outwards 
from the disc plane. This means that in the jet approaching the observer B^ = —fl F rB z /c, whilst in the jet receding 
from the observer B^ = +fl F rB z /c. Disturbances propagating along the approaching jet should be brighter than 
those in the counter-jet due to the Doppler boosting effect. Therefore, we will observe knots in parts of jets having 
B^ = —Q F rB z /c more often than in those parts where B^ = -\-Q F rB z /c. This can explain the fact that the number 
of sources having II > is somewhat larger than the number of sources having II < despite the fact that for most 
of the values of the parameters calculations give a longitudinal polarization II < 0. We see from Fig. [S] that our model 
curves fit the data for BL Lac objects better than those for quasars. On the other hand, the polarization measurements 
for the quasars are quite different at different frequencies and such differences are larger than the variations of the 
polarization measurements for the BL Lac objects. In the frame of our model, jets in quasars tend to have higher 
values of VL F than do jets in BL Lac objects. It is remarkable that in quasars knots in jets also have faster apparent 
motions than in BL Lac objects. 



4. Summary 

We calculate the linear polarization of the synchrotron radiation produced by fast particles in non-uniformly moving 
and differentially rotating relativistic flow such as jet flows in AGNs. The magnetic field in the jet has axial and 
toroidal components. The distribution of the fast particles is taken to be isotropic in the reference frame comoving 
with the flow. We perform integration of Stokes parameters /, U, Q along the line of sight passing through the jet 
with all relativistic effects fully taken into account. We consider two possibilities for the spatial density of the fast 
particles: 1) ultra-relativistic particles distributed uniformly across the jet, 2) ultra-relativistic particles concentrated 
close to the Alfven resonance surface. The latter distribution follows from our acceleration theory and we compare it 
with the former which might be produced by other acceleration processes. We assume a uniform axial magnetic field in 
the jet, so the equilibrium configuration of the jet depends only on the radial profile of the angular rotational velocity 
£l F (r) of the magnetic field lines. We take the profile of the angular rotational velocity to be a plausible monotonic 



function of radius, O 




and investigate the dependence of polarization on the magnitude of the 



10 



V.I. Pariev et al.: Relativistic parsec-scale jets: II 



angular velocity f2. In our computations we also specify a component of the velocity parallel to the magnetic field lines 
such that the kinetic energy associated with the flow is minimized. The actual velocity is determined by the process 
of jet formation, which is poorly understood now; still such minimal velocity is physically favoured (Frolov & Novikov 
1998). We calculate the degree of linear polarization and orientation of the electric vector with respect to the jet axis 
for the emission integrated over the jet cross section, taking into account the relativistic bulk spiral flow of plasma in 
the jet. This flow leads to the presence of an electric field comparable in magnitude with the magnetic field, causes a 
Doppler boost of the intensity of the synchrotron radiation, and a swing of the polarization position angle. We want 
to stress that the last effect, which was first pointed out in 1979 by Blandford & Konigl for the simpler case of a 
uniform magnetic field parallel to the jet axis, substantially modifies the resulting degree of linear polarization from 
the whole jet. However, because of the axial symmetry of the problem, only two orientations of polarization position 
angles are possible: parallel or perpendicular with respect to the jet axis. Observations of polarization of VLBI scale 
jets in BL Lac objects on centimetre wavelengths (Gabuzda et al., 2000) show a bimodal distribution of the position 
angles: most of the knots in BL Lac objects have polarization either parallel to the jet axis or perpendicular to it. 
Similar observations of quasars (Lister 2001) on millimetre wavelengths indicate a unimodal distribution of the position 
angles: most of the knots in the observed quasars have the electric vector of polarized emission oriented along the jet 
axis. In our present calculations the polarization does not depend on frequency since we have assumed a power-law 
distribution over the whole range of energies of the particles. Different slopes of the energy spectrum of the particles at 
different energies can result in a dependence of the degree of polarization and of the position angle on the frequency. 
This may account for the differences in polarizations at different frequencies observed in quasars. At present, we can 
only suggest that jets in quasars tend to have larger f2 F and a larger ratio of toroidal to poloidal magnetic field than 
do jets in BL Lac objects. This can account for the electric field vector in quasars being predominantly parallel to the 
jet axis. 

We also calculated the velocities of crests of standing eigenmodes and compared them with the observed proper 
motions of bright VLBI knots in jets. We were able to match the observational data only for the case when the 
emitting particles are concentrated close to the Alfven resonance surface and the magnetic field lines are twisted in 
the direction opposite to the direction of rotation f2 F . The last fact can be understood by a selection effect due to the 
Doppler boosting which makes jets approaching the observer more favourable to observe than receding jets. Comparing 
calculations with observations we can estimate the angular rotational velocity of the jets. The trend is that tt F in 
BL Lac objects is intrinsically smaller than in quasars. More detailed calculations considering different possibilities 
for the £l F (r) profile are postponed to another work. There are some discrepancies between our calculations and the 
observation points (see Fig. EJ. It is possible that the radial profile of fi F (r) needs to be tuned to provide better 
agreement. At present it is difficult to suggest a unique test, which would make it possible to distinguish between 
electromagnetic and shock mechanisms for the origin of knots in jets. Future high-resolution and high-dynamic-range 
observations could provide detailed images of the spatial structure of knots and surrounding jets, and provide the 
answer. 

Acknowledgements. The authors are grateful to D.C. Gabuzda for compiling observational data on polarization of jets and 
helpful discussions of different observational issues. This work was done under the partial support of the Russian Foundation 
for Fundamental Research (grant number 02-02-16762). Support from DOE grant DE-FG02-00ER54600 is acknowledged. 



V.I. Pariev et al.: Relativistic parsec-scale jets: II 



11 



References 

Attridge J.M., Roberts D.H., & Wardle J.F.C., 1999, ApJ, 518, L87 

Begelman M.C., Blandford R.D., & Rees M.J., 1984, Rev. Mod. Phys., 56, 255 

Beresnyak A.R., Istomin Ya.N., & Pariev V.I., 2003, A&A, accepted 

Blackman E.G., 1996, ApJ, 456, L87 

Blandford R.D., & Konigl A., 1979, ApJ, 232, 34 

Blandford R.D., & Eichlcr D., 1987, Physics Reports, 154, 1 

Cocke W.J., & Holm D.A., 1972, Nature Phys. Sci., 240, 161 

Frolov V.P., & Novikov I.D., 1998, Black hole physics: basic concepts and new developments. Kluwer Academic, 
Dordrecht 

Gabuzda D.C., & Cawthorne T.V., 1993, in Davis R.J., Booth R.S., eds, Subarcsecond Radio Astronomy. Cambridge 
Univ. Press, Cambridge, p. 211 

Gabuzda D.C., Mullan CM., Cawthorne T.V., Wardle J.F.C., & Roberts D.H., 1994, ApJ, 435, 140 
Gabuzda D.C., 1999a, New Astronomy Reviews, 43, 691 

Gabuzda D.C., 1999b, in Ostrowski M., Schlickciscr R., eds, Plasma Turbulence and Energetic Particles in Astrophysics. 
Uniwersytet Jagiellonski, Cracow, 1999, p. 301 

Gabuzda D.C., 2000, in Conway J.E., Polatidis A.G., Booth R.S., Pihlstrom Y.M., eds, EVN Symposium 2000, 
Proceedings of the 5th European VLBI Network Symposium. Onsala Space Observatory, Onsala, p. 53 
Gabuzda D.C., Pushkarev A.B., & Cawthorne T.V., 2000, MNRAS, 319, 1109 
Gabuzda D.C., Pushkarev A.B., & Garnich N.N., 2001, MNRAS, 327, 1 

Ginzburg V.L., 1989, Applications of Electrodynamics in Theoretical Physics and Astrophysics. Gordon and Breach 
Science Publishers, New York, p. 65 

Hirotani K., Iguchi S., Kimura M., & Wajima K., 1999, PASJ, 51, 263 
Hirotani K., Iguchi S., Kimura M., & Wajima K., 2000, ApJ, 545, 100 
Hutchison J.M., Cawthorne T.V., & Gabuzda D.C., 2001, MNRAS, 321, 525 
Istomin Ya.N., & Pariev V.I., 1994, MNRAS, 267, 629 
Istomin Ya.N., & Pariev V.I., 1996, MNRAS, 281, 1 
Konigl A., & Choudhuri A.R., 1985, ApJ, 289, 188 
Laing K.A., 1981, ApJ, 248, 87 

Leppanen K.J., Zcnsus J. A., & Diamond P.J., 1995, AJ, 110, 2479 
Lister M.L., Marscher A.P., & Gear W.K., 1998, ApJ, 504, 702 
Lister M.L., 2001, ApJ, 562, 208 

Lister M.L., Tingay S.J., & Preston R.A., 2001, ApJ, 554, 964 
Manolakou K., Anastasiadis A., & Vlahos L., 1999, A&A, 345, 653 
Morrison P., & Sadun A., 1992, MNRAS, 254, 488 
Pelletier G., & Pudritz R., 1992, ApJ, 394, 117 

Reynolds C.S., Fabian A.C., Celotti A., & Rees M.J., 1996, MNRAS, 283, 873. 

Romanova M.M., & Lovelace R.V.E., 1992, A&A, 262, 26 

Roser H.-J., & Meisenheimer K., 1991, A&A, 252, 458 

Sazonov, V.N., 1969, AZh, 46, 502 

Sikora M., & Madejski G., 2000, ApJ, 534, 109 

Scheuer P.A.G., 1984, Adv. Space Res., 4, 337 

Wardle J.F.C., Homan D.C., Ojha R., & Roberts D.H., 1998, Nature, 395, 457 
Zhelezniakov V.V., 1996, Radiation in astrophysical plasmas. Kluwer, Dordrecht 



12 



V.I. Pariev et al.: Relativistic parsec-scale jets: II 




Fig. 1. The scheme illustrating the geometry of the jet and of the polarization observations. The jet boundary for 
r = R and three magnetic tubes for r = i?/4, 2i?/3 and 9i?/10 are shown for uniform B z and angular velocity 
profile ljl8|l with fl = 10. The magnetic field lines are spiralling on a magnetic tube. Since tt F (R) ~ 0, the total 
current through the jet is equal to zero and the magnetic field is purely poloidal both at the boundary and at the axis 
of symmetry. The density of the poloidal current j z is negative when r < R/V2 and positive when R > r > R/y2. 
The electric field E induced by jet rotation is radial. The plasma velocity v along the magnetic tube consists of two 
components: rotation with angular velocity £l F (r), and motion along the magnetic field lines with a speed V| = KH. 
A line of sight OQ passes through the jet at a closest distance CP = h from the jet axis. The line CA is parallel to 
the plane containing the observer and the jet axis. OXY is a rectangular coordinate system on the plane of the sky of 
the observer, OX is in the plane containing the observer and the jet axis. The line OB is parallel to the jet axis, and 
the wave vector of a photon k makes an angle 9 with the line OB and the jet axis. Vector s is a unit vector in the 
direction of the electric vector of the polarized component. Vector s is in the plane OXY and makes an angle \ with 
OX. 
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Fig. 2. Dependences of linear polarization II on angular rotational velocity (left panel) and on visual velocity of 
knots v V [ S (right panel) for fundamental axisymmetric mode m = 0. II > when the electric field vector is oriented 
parallel to the jet axis, II < when the electric field vector is oriented perpendicular to the jet axis. v v - m > means 
that the observed motion of knots is direct, d v ; s < means reversal motion. Crosses are for the BL Lac objects and 
squares are for the quasars listed in Table 1. 



-0.5 





Fig. 3. The same as in Fig. |21 but emitting particles are concentrated close to the Alfven resonance surface, m = 1 
fundamental mode. The curves are plotted only for 0.7 < < 4.9, when the Alfven resonance exists for the m = 1, 
k = fcmin mode. 
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Fig. 4. Dependence of the polarization II on (left panel), on v v [ s (middle panel) for £ = 2.1 and for £ = 3 (right panel). 
The emitting particles are concentrated close to the Alfven resonant surface. Here B$ — +fl F rB z /c, fundamental 
eigenmode with m = 1. Notations are the same as in Fig. [3 In the left panel solid curves are for ( — 2.1, dashed 
curves are for £ = 3. 




Fig. 5. The dependence of polarization II on (left panel), on v V [ S for £ = 2.1 (middle panel) and for £ = 3 (right 
panel). Emitting particles are concentrated close to the Alfven resonant surface. = —il F rB z /c, fundamental 
eigenmode with m — 1. Notations are the same as in Fig. In the left panel solid curves are for ( = 2.1, dashed 
curves are for ( = 3. 



